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ABSTRACT 

The flow field inside a whirling annular seal has been measured using a 3-D Laser Doppler Anemometer 
(LDA) system. The seal investigated has a clearance ot 1 .27 mm, a length of 37.3 mm and is mounted on a drive 
shaft with a 50% eccentricity ratio. This results in the rotor whirling at the same speed as the shaft rotation 
(whirl ratio = 1.0). The seal is operated at a Reynolds number of 12,000 and a Taylor number of 6,300 (3,600 
rpm). The 3-D LDA system is equipped with a rotary encoding system which is used to produce phase averaged 
measurements of the entire mean velocity vector field and Reynolds stress tensor field from 0.13 mm upstream 
to 0.13 mm downstream of the seal. The mean velocity field reveals a highly three dimensional flow field with 
large radial velocities near the inlet of the seal as well as a recirculation zone on the rotor surface. The location 
ot maximum mean axial velocity migrates from the pressure side of the rotor at the inlet to the suction side at 
the exit. Turbulence production is a maximum near the seal inlet as indicated by the rapid increase of the 
turbulence kinetic energy (/c). However, turbulence production and dissipation attain equilibrium fairly quickly 
with k remaining relatively constant over the last half of the seal. 

NOMENCLATURE 

A = leakage area, tt-Dc, m : 

c = nominal clearance between the rotor and stator = R2-R1, m 
D = stator diameter, m 
L = rotor length, m 

0 = leakage rate, m 3 /s 

Re = Reynolds number = 2p Uc/'n 

R1 = rotor radius, m 

R2 = stator radius = D/2. m 
Ta = Taylor number = [/?W sh c/|j][2c/D] " 

U = average leakage velocity = Q/A, 3.71 m s 
fT = axial mean velocity, ms 
tT r = radial mean velocity, m/s 
Uq = azimuthal mean velocity, ms 
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axial time averaged Reynolds normal stress, m 2 /s 2 
radial time averaged Reynolds normal stress, m 2 /s 2 
azimuthal time averaged Reynolds normal stress, m 2 /s 2 
azimuthal velocity of the rotor surface, 30.9 m/s 
axial distance from seal inlet, m 

azimuthai angle measured from the minimum clearance in the direction of shaft rotation, clockwise 
on the figures 

turbulence kinetic energy = V 2 (u ; u ; +tT, u ; + uju 9 ), m 2 /s 2 
absolute viscosity of the fluid (water), Ns/m 2 
fluid density, Kg/nr 

fluctuating value = instantaneous value - mean value 
mean value (time averaged) 


INTRODUCTION 


Hie influence of seals upon the rotordynamic stability of turbomachines has become an important concern 
in the last 20 years as turbomachines have become lighter and have been required to operate at increasingly 
higher speeds. This was clearly demonstrated on the Space Shuttle Main Engines (SSME) where the seals became 
a critical element in the stable operation of the turbomachine. Accurately predicting the leakage rate through 
a seal is an extremely difficult task due to the great number of geometric variables as well as harsh operating 
conditions and exotic fluids. Such complexity pushed the Limits of empirical and bulk flow models. Now with 
the additional requirement of knowing the forces generated by the seal and how they effect the rotordynamic 
stability, a new generation of analyses is needed. 

To address this challenge, Tam, Przekwas, Muszynska, Hendricks, Braun, and Mullenfl] undertook a 
numerical and analytical study of the fluid dynamic forces in seals and bearings. They solved the Reynolds 
averaged Navier-Stokes equations for the turbulent, swirling flow field inside annular seals using the PHOENICS- 
84 numerical code. In their st u dy[ 1 ], the effects ol preswirl, whirl ratio, overall pressure drop, shaft speed, fluid 
injection, and eccentricity ratio upon the flow field were examined and an evaluation of the resultant fluid 
dynamic forces generated by the seal was performed. These calculated forces and overall leakage rates were then 
compared with measurements performed by Childs[2,3]. Overall, good agreement was obtained. 

OBJECTIVES 


There is one obstacle standing in the way of further advancement of this type of computational analysis - 
the turbulence model. The statement of Tam, et.al.fl], "until fluctuation details for high shear flows at low- 
clearances become available the simplest of the turbulence models was selected for the present calculations," 
indicates a need for more detailed basic information. In order to adequately evaluate the effectiveness of the 
turbulence model, details of the flow field (velocity field and the Reynolds stress tensor) are required. It is the 
objective of an ongoing research project at the Turbomachinery Laboratory, funded by NASA Lewis, to obtain 
this required data. Presented in this paper is a portion of the results. In particular, the distributions for the three 
mean velocity components and the turbulence kinetic energy level are presented for one axial location 
immediately upstream of the seal inlet and five axial positions along the length of the seal. Additional 
information, including the Reynolds stress tensor, was obtained at various other axial locations; space limitations, 
how'ever, prohibit the presentation of all this information in one paper. The data are available on MS-DOS disks 
and a more complete presentation is available in Mr. H.D. Thames’ thesis[4], 

FACILITIES AND INSTRUMENTATION 

Hie facility used in this investigation has been previously described in detail by Morrison[5,6]. Briefly, 
the test section (Figure 1) consists 50.4 mm diameter overhung shaft directly coupled to a 37 kW electric 


induction motor. The electric motor is driven by a power supply capable of varying speeds from 200 to 5,300 rpm. 
The seal itself is composed of an acrylic rotor (0. 164 m O.D.) and a stainless steel stator with a nominal clearance 
between the rotor and stator that is 1.27 mm. The rotor is optically coated by a vacuum-deposition technique 
to reduce reflected light intensity and is mounted on the stainless steel shaft using a brass bushing, the outer 
circumference of which is 0.63 mm eccentric with reference to the axis of the shaft. This results in the rotor being 
eccentric causing whirl in a circular orbit at the same speed as the shaft (whirl ratio 1.0). The seal stator is 
constructed of stainless steel with a smooth finish. To permit access for the LDA system, a flat optical window 
is installed in a narrow slit along the axis of the seal. A flat window is used to avoid lens effects. The deviation 
from roundness resulting from the window is only 1.5% of the clearance or 0.027% of the rotor radius. The 
annular seal is a constant radius design which is 37.3 mm long. 

A 3.75 m 3 storage tank supplies water at approximately 300 K to a centrifugal pump. The pump flow rate 
of 0.00243 m'7s is set by a throttling valve and metered by a turbine flow meter. After metering, the water is 
introduced into the seal rig and water exiting the seal is returned to the storage tank. Expanse! 462 WU, 6 pm 
diameter plastic sphere with a specific gravity of 1.23, is added to the water to provide light scattering seed for 
the laser Doppler anemometer system. This particle will follow instantaneous flow fluctuations up to 44,000 Hz. 

Morrison, et al [5,6] have previously discussed in detail the 3-D LDA system used and have presented 
an error analysis for the system. The LDA system is a three color, six beam arrangement with 8.5X beam 
expanders. Two colors measure orthogonal velocity components and the optical axis for the third color is inclined 
at 30° to the others. This results in a measurement volume 0.025 X 0.025 X 0.10 mm in size. Bragg cells are 
installed on one beam of each color pair so that flow reversals can be measured. This LDA system is thus 
capable of measuring the instantaneous 3-D velocity vector at the measurement volume. The data acquisition 
systems consists of a rotating machinery resolver and three counters which supply the angle of the rotor as well 
as the three non-orthogonal velocity components to an 80386 computer system. A gate on the master interface 
assures that all three counters measure a velocity component within a coincidence window of 10 ps. The data 
is then post processed to calculate the mean velocity, Reynolds stress tensor, and other quantities with the data 
divided into 20 different angular positions of the rotor, each with a window of 18° (one window had inadequate 
data for analysis). A minimum of 90,000 velocity realizations sorted into the different angular positions of the 
rotor were recorded at each X and r location. This resulted in a minimum of 1000 and a maximum of 
approximately 8000 samples at each reported location for a given phase window. Figure 2 presents the locations 
on the r-0 plane where the measurements were made. 

Recent experience (Wiedner[7] and Morrison[8]) with this system has shown that there is some velocity 
bias present in the data which is best corrected using the McLaughlin-Tiederman[9] velocity bias correction 
technique. Uncertainties are estimated[5] to be ±1% on and £7^, ±2% on and 20% on tc. 

RESULTS 

Isovels were calculated from the data for the axial (£7^), radial (£7^), and azimuthal (£7^) mean velocity 
components at X/L = -0.0018, 0.00, 0.21, 0.50, 0.77, and 1.00. Plots of these are presented in Figures 3-5. The 
rotor is orbiting in the clockwise direction (increasing 8). 

Hie seal influences the mean velocity field immediately upstream (X/L = -0.018) of the inlet in two ways. 
First, because the large end of the rotor is exposed to the supply plenum, there is an induced preswirl which varies 
from 6% of the rotor surface speed (0.06W lh , 0.50U) near the stator wall to 28% near the rotor. Secondly, the 
flow accelerates into the seal with the maximum axial velocity (1.3U) occurring near 8 = 105° rather than at the 
maximum clearance location (8 = 180°). With the rotor spinning clockwise, it appears that fluid approaching the 
inlet at the region of maximum clearance maintains the axial momentum necessary to enter the seal but must 
accelerate as it is rotated into the smaller gap. This hypothesis is further substantiated in the analysis of the axial 
velocity contours at X/L = 0.00 where the fluid is semi-confined yet the precession of the rotor causes the axial 
velocity to increase even more (2.4U) as the clearance decreases ( 8 ~ 75°). It is interesting that the fluid not only 
continues flowing downstream as the clearance decreases but actually experiencesthis substantial increase in local 
axial velocity. This is consistent with continuity since the azimuthal velocity remains small compared to the axial 
velocity at this azimuthal location. In this region, large outward radial velocities (0.9U) are caused by the rotor 



moving toward the stator. In fact, the (T r lT 9 velocity vector in this plane is almost entirely in the outward radial 
direction for the rotor phase 30° <. 8 <_ 180°. This combined with the steep dUjdi near the rotor indicate the 
presence or onset of a flow separation on the rotor (as evidenced by the close spacing of the contours). 

On the opposite side of the rotor (the phase where the clearance is increasing, 8 > 180°), the mean axial 
velocity decreases as 8 increases until there is a negative axial velocity component (flow from the seal back into 
the plenum) from 8 = 270° to 350°. The minimum velocity is -0.4U at X/L = -0.018 and -0.3U at X/L = 0.00. 
This side of the rotor (180° < 8 < 360°) is typically called the suction side. Therefore, one might expect this 
region where the clearance is increasing to provide less resistance to flow and to "suck in" fluid which would result 
in a positive C x . However, in this case, a strong azimuthal velocity component, is bringing fluid into the 
widening clearance at a faster rate than the expansion can support causing a pressure increase. This results in 
the fluid actually 'squirting out" of the seal back into the inlet plenum. 

Progressing from XT = 0.00 to 0.21, the region of maximum moves counter-clockwise to 8 = 90° and 
outward radially while decreasing in maximum value from 2.4U to L5U. The region occupied by the high speed 
flow is more uniform with smaller dVjd r gradients near the rotor indicating the reattachment of the flow to the 
rotor. U x is positive across the entire cross-section with a minimum value of 0.2U, which indicates that a 
stagnation zone exists somewhere between X/L = 0.00 and 0.21. The radial velocities have become less than 10*2 
of U compared to 80% at X/L = 0.00 indicating the flow is becoming more parallel to the stator. The azimuthal 
velocity distribution is much more uniform around the rotor being almost symmetric about 8 = 0°. 

Bv the center of the seal, X/L = 0.50, £7] has been redistributed with the decreasing maximum velocity 
shifting clockwise. The flow is becoming established with the mean axial velocity being very nearly uniform over 
the entire seal and the radial velocity decreasing even further in amplitude to 0.04U. The azimuthal velocity 
distribution remains very much the same as at XT = 0.21 except that there is a slight overall increase in 
magnitude. 

The axial velocity component, IT V continues to show development by actually rotating locations of 
maximum value (1.5U) to 6 ~ 315°, where the backflow existed in the inlet, and the minimum £T (0.4U) occurring 
where the maximum was present at the inlet. The radial velocities have increased slightly with a longer cross 
section of the flow having + 0.04U. In the transition from X/L = 0.21 to 0.77 the regions of positive and negative 
radial velocity have exchanged places. The azimuthal velocity has increased slightly since X/L = 0.50 attaining 
approximately 35% of the shaft speed for 8 < 90° and 25% for 8 > 200°. The azimuthal velocity distribution is 
the driving mechanism for the change in location of maximum mean axial velocity. From XT = 0.50 to the end 
of the seal, U 9 is consistently larger for 8 < 180°. This causes a migration of the fluid with the maximum (T x from 
8 ~ 75° at the inlet to 8 ~ 315° at XT = 0.77. This effect continues until the exit of the seal where U[ attains 
the maximum value of 1.9U and the minimum value of 0.3U. In the exit plane the radial velocities begin to rise 
in value as the flow begins to rearrange for the step change in area at the exit. 

One of the variables calculated in most turbulence models is the turbulence kinetic energy 
(zc = V:{ u s u t + u l u , + u 9 u ' 9 }). It is a good measure of the level of turbulence present in the seal. As mentioned 
earlier, we have measured the entire Reynolds stress tensor which is available in Thames[4]. The contours of k 
presented in Figure 6 show that immediately upstream of the seal, the turbulence level is large with a maximum 
k = 0.4U ; . This corresponds to a turbulence intensity which is approximately50% of the average leakage velocity 
(U). This occurs near the rotor at 8 = 270° where is zero. Over the majority of the field, tc ranges from 0.2U* 
to 0.3U : (35% to 45% turbulence intensity), however, once the flow enters the seal, the turbulence kinetic energy 
increases rapidly in the region near the rotor where the mean flow appears to separate, and hence k exceeds 
0.8U : . This is consistent with the very large radial gradient of the mean axial velocity (dUjd r) which contributes 
to the production of turbulence. As the flow progresses further downstream the turbulence kinetic energy levels 
decrease to approximately 0.3U 2 at XT = 0.21 and 0.2U : for 0.50 <. X/L <. 1.00. This relatively gradual change 
is expected since the spatial gradients of the mean velocity are significantly smaller for xT >_ 0.21 hence the 
mechanism for turbulence production is reduced and turbulence dissipation begins to balance the turbulence 
production. 


CONCLUSIONS 


The flow field inside a whirling annular seal has been measured using a 3-D LDA system. The mean 
velocity field reveals a highly three dimensional flow field with large radial velocities near the inlet of the seal and 
a recirculation zone on the rotor surface. Near the inlet the mean axial velocity is very large on the pressure side 
of the rotor but by the downstream half of the seal, the region of high mean axial velocity has migrated to the 
suction side. Near the inlet, on the suction side of the rotor, there is a back flow into the plenum which is caused 
bv the large azimuthal velocity moving large quantities of fluid into the region which then ’’spill over’’ into the 
plenum. The radial velocities steadily decrease in value as the flow progresses through the seal until the exit 
plane is reached at which point they increase abruptly, however, the magnitude is still significantly lower than at 
the inlet. The overall azimuthal velocity steadily increases from approximately O.IZW^ at the inlet to 0.30W sh at 
the exit. The flow over the seal inlet produces very large levels of turbulence kinetic energy. However, 
turbulence production and dissipation attain equilibrium fairly quickly and k remains relative constant over the 
last half of the seal. 
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